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ABS Gane e > 9 gre 


A control system was designed to attenuate vertical 
accelerations for the XR-3 captured air bubble type surface 
effect ship using linear regulator techniques applied to the 
Simplified nonlinear equations of motion. A pressure lift- 
only model was used to represent the craft vertical heave 
motion and was linearized around the steady state operating 
point. Model validation was obtained through analysis of 
the frequency spectrum. State variable feedback was used to 
determine a set of optimal control gains that would reduce 
the magnitude of the heave acceleration during operation 


under simulated sea input conditions, 
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I, INTRODUCTION 


An important design consideration for the high speed 
captured air bubble type surface effect ship is the vertical 
acceleration associated with operation in a random ocean 
environment [Ref, 1]. The vertical motion, heave, is caused 
by waves changing the volume of the plenum producing pres- 
Sure variations that result in vertical craft accelerations. 

Studies of conventional hull ships show that the strongest 
vertical plane accelerations are experienced at encounter 
frequencies close to the ship's natural frequency. Reference 
2 states that "The most severe and objectionable motions (wet 
decks, slamming, high accelerations) result from the responses 
to the components that are near synchronism with the ship's 
natural pitching and heaving periods". Encounter frequency 
is a function of wave frequency, craft speed and direction. 
Since most of the weight of the surface effect ship 1s sup- 
ported by a bubble of air, large accelerations can occur that 
affect habitability and sea keeping ability while operating 
in high seas. 

Methods of attenuating vertical motion have been the 
subject of past studies. A conventional design by Rohr Marine 
Inc. [Ref. 3] utilized accelerometer output to monitor verti- 
cal motion variations. These signals are passed through 


compensation circuits, filtered and applied to actuators 








which control air flow rates into the plenum. State variable 
techniques are applied to describe the system and control 
functions. The physical constraints of the actuators are 
considered with regard to their inherent nonlinearities and 
limitations, The sea input is modeled using a wave height 
power spectral density procedure, 

Through the use of a towing tank model and data gathered 
from the 100B test craft, an investigation was conducted to 
measure the characteristic response of the SES while operating 
in an actual sea state environment [Ref. 4]. The model was 
subjected to simulated sea state inputs of various severity in 
an effort to establish the linear characteristics of the 
vertical motion and determine the damped natural frequency of 
response. A series of curves showing response amplitude for 
both C.G. accelerations and pitch angle versus encounter fre- 
quency were developed for a range of speed and sea state 
operation. 

In Ref. 1, the frequency response approach was use¢d in an 
analysis which examined the effects of changing plenum air 
flow rates on center of gravity accelerations using the six 
degree of freedom simulation program of the XR-3 [Ref. 5]. 
iieeasseudy deyeloped a linear simplified model in close agree- 
ment with the 6 D.O.F. model and in addition used signifi- 
cantly less computer time, 

Tn Ref. 6, Boggio conducted a study on the effect of the 
use of a flexible membrane in the plenum chamber itself. 


Using the loads and motions program for the XR-3, the study 








reduced computation time through the inclusion of pressure 
rate and volume rate equations. Membrane insertion served to 
effectively reduce the magnitude of pressure variations in 
the plenum thereby attenuating vertical motions. 

Grant [Ref, 7] used simple models for the control devices 
applied to the loads and motions program of the 100B SES 
craft. From monitored plenum pressure variations and heave 
accelerations, control signals were developed by filter de- 
Sign methods, to vary the fan speed and vent louver openings 
in order to reduce the vertical accelerations and pitch motion. 

This thesis applies the linear regulator design technique 
to determine a set of optimal state feedback gains to control 
the air flow rate. The equations of motion used in the design 
were derived from a simplified pressure lift only model lin- 
earized around the steady state operating point. Design re- 
sults were evaluated using both the linear and nonlinear 


models. 





Pe RESEND ALL ON OP eres ofoleM 


Two major forces which support the weight of the surface 
effect ship are the buoyancy and plenum pressure forces. 
The buoyancy forces (Fi) are exerted on the craft by dis- 
placement of water by the underwater volume of the rigid 
Sidewalls. The plenum piasegunce lift force (F,) is exerted 
on the craft by the captured cushion of air. Flexible bow 
and stern seals along with the rigid sidewalls and lift fans 
form the captured air bubble. Relative contributions to 


craft vertical plane motion by Py and F, depend upon the 


Ie 


encounter frequency (wi). At low frequencies, F, responds 


L 
to the rise and fall of the sea. At high W > F, attenuates 
because incremental volume changes of the immersed hull are 
reduced. The Po force responds inversly to changes in plenum 
volume (V,) according to the adiabatic law of pressure 
variations within the plenum chamber [Ref. 1]. Volumetric 
changes in the plenum depend upon the sea state and encounter 
frequency. As is shown in Ref. 1, the effect of the volume 
changes produces an increase in the Fo force as the encounter 
frequency is increased, 

In the system developed in this study, the buoyancy terms 
have been omitted in an effort to simplify the system further. 


This simplified model is referred to as the "pressure lift 


enmy esystem. Transfer functions of acceleration, volumetric 


0 





plenum air mass to input waveheight and system state equa- 
tions were developed from the linearized form of the pres- 
ere slaiit Only system. 
The basic assumptions which form the basis for the pres- 
sure lift only model are presented below: 
1. The center of pressure ee is directly under the 
Comecm or gravity (C.G,). 
2. The sidewalls are of uniform cross section and are 
symmetrical about the C.G. 
3. The effect of pitch moments on the vertical plane 
forces are neglected. 
4. The only variable lift forces are those forces 
caused by changes in plenum pressure (P.). 
5. The craft 1s at a constant speed. 
6. The air leakage rate from the rear seal is constant. 
7. The craft weight 1S supported by two component 
forces in the following approximate proportions; 
Sidewalls: 20% (fixed force) 
Captured Air Bubble: 80% (variable force) 
8. Encounter frequencies of interest lie within the 


range of 17<w, <100 radians/sec. 


A, EQUATIONS OF MOTION 

The development of the nonlinear simplified model for 
vertical motion is presented in Ref, 1 and the block diagram 
appears in Appendix A. The development of the pressure lift 
only nonlinear model follows in a similar manner. The equa- 
tions of motion differ in that the F, force is considered 


L 


Gonstant. 


yeah 





The application of the given assumptions to the Newtonian 


Laws for vertical motion of a rigid body system leads to the 


following equations; 


1) 


2) 


ap 


4 ) 


5) 


6 ) 


Orifice Leakage Rate 


q oe Cage 2 en aoe 

out Tne ale Os sec 
Fan Map Input Rate 

z ft 

eRe = nlQ, (a) - KP pd se 


Absolute Plenum Pressure 





Db lbs 
res Pe ) 
b a Hess et? 
Plenum Air Flow Rate 
lbs 


7, = (W = ‘p = =e ft? 
ee M M sec 


Heave Velocity 


a a3 ipa 
Zs fi alit aaa 


eZ 


CZ. 


Cee 


C2 


CZ 


(2. 


C2 


oe 


Pee 


te.) 


4) 


oS) 





7) Heave Displacement 

n= f Z at ft (2.4.7) 
8) Plenum Pressure Lift Force 

ego A. P lbs C27 fc) 
eesuoyancy Lift Force 


Bay = [2A 1 J ibs C27. 9) 


it d?p 


where 1, is steady state draft. 


d 


10) Plenum Gage Pressure 


Po = (P, - P_) — (2.A.10) 
a 
eye pratt 
Ae ee (2.A.11) 
d S 


The system parameters and constants are listed as follows: 


1) Adiabatic Process Coefficient 
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Zz) 


3) 


4 ) 


3) 


6 ) 


> 


8 ) 


Leakage Area 


Air Density 


ee mUOOZoTS 


Atmospheric Pressure 


ees eS 
a 


Plenum Area 

A. = 200 et? 
b 

Empty Plenum Volume 

V = 383 Ft? 
n 

Craft Weight 

We= 6007 Eos 


14 





9) Area of Each Immersed Sidewall 
pee 7 Of US 18 75 Tic 


mpm pratt Initial Condition 


L Co) = ,9 2 

11) Density of Water 
9 =1.99 

Meee Gravity Acceleration Constant 
G = 32.174 ft/sec“ 


13) Vertical Center of Gravity Location Above Keel 


bane 
Fee 2.5 feats 


ieee Number of Fans 


15) Steady State Fan Output 


Q; Co) = 35 ft°/sec 


ses 
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PeeeolNEAR MODEL 
Prior to defining the system states, the non linear 
equations were linearized around their steady state operating 
conditions by application of a Taylor series expansion. The 
details of the linearization process is shown in Appendix A. 
The resulting simplified pressure only signal flow graph 
for the linearized system is shown in Figure 2. 


1. Transfer Function Development 


The system graphical determinent is; 


ee (Ly a L.) + Ly L, C2 ea. 


hee lea) (L.-+ L.) ee.) 


where 
— _-Ll_-1 
= S Ss (-kg)(-k,) (k,) ¢-k,) 
5G, 


= Fa N ees) 
S kyk 3k, ky ( ) 


ee 
L, = 8s (ko) Ck -k,) 


1 ena) 
-S kk. ko , 


The transmission gain from the control input to the accel- 


eration output for calm water operation is; 
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G@ odnst yj 


Oa A 
SUOT]TDUOOD = 
SUT} eJedO Ja eM WTPRO UOJ es aed = (ie 
XNTJUT ute sntd 33d aseyeoelt Gee Zo 
AEB doy wUSTOTJZe0o |yz st Py- :aLON Teele || Za 





NOLLVasd TIOOV 
LAG LNO 


€ > HdVu¥9 MOTTA 
TVNOIS ATINO L421 
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ILS, 





9-72 (248.5) 


where 
P_ = -k.k, k-k..s7- (D> 
al Se eae le ene 
and 
Ay = 1 Coe) 
and 
a, = 0 C2 Bie8) 
Pest bubing equations (2.B.2), (2.B.6), (2.B.7) into 
equation (2.B.5) yields; 
epee 3 
ae 2. 3°47 "10 u 
a, = 0 ee a ns een C223 9.) 
1-[-s ok3kykg)+C-s k,k,k) J 
Redwerme further; 
-k.k, kk f) 
a, = te (243210) 
S tk,k. ko Stk kK, Ke, 


The transmission gain from the sinusoidal sea generator input 
to the acceleration output with no control signal is; 


Raa 


ee ce 
a, = Se ‘Oy (2s a) 


20 





where: 


P, = -k)k,k,k, Caen) 
and 

Ay ak C2 opal 3.) 
and 

oF = 0 CZ. bY) 

Appropriate substitution into equation (2.B,.11) 
yields 

-k,kyk,k, WW 

Ze -2 —1 

1-[(s “kpk,k,kg)+(-s “k3k,k,) J 


Reducing further, the transfer function is 


: 2 
2, | Sep ESe sh Os 
ow 24 se a a. — 
S tkikpk7Stkok3k, Kg 


For the condition that control input is applied, the 


acceleration output has the following transform; 


i ak kkk, 984, 
- P= nee cast 
Zi 
S tk kk Stk, kkk, 


2 
kK, kk,k,s° Ww 


_ +, 


oe C273 4169 


2 
S tk, kp k7stk) kk, Kg 





Zale 





Following the same procedure for the applied control 


input, the volumetric plenum air mass output has the 


following transfer function; 


3 1 
5 -1 
P = Kis 
al 
> a aI Solel 
k. .stk,k.k, k ao 
ny = 10 pas eat 4 
™ 2 u 
S #k, kk stk k,k ky 


Pamerrom (2.8,15) is of the form; 


OUTPUT _ 
new yee 


Ero 


+bstc 


and with no control action, becomes; 


y 
_y(s) = ae 
ow (Cs ) eee 
where 
y(s) = 22 
and 
By 6s) b 





Ws. Bee 


Ze 


(Coe 


wd.) 
G7 as) 
alee) 
Cao. 20) 
G2 Bc 


2A); 


Coe B22} 


(2258423) 


24) 


223.2) 





pileimer tne Coefficients consist of the following parameters; 


a -2-k) k,k3k, C7825) 
b= kjk. ko 725.27) 
e = kjk 3K Ke C223, 29) 
d =~k, KK 3k, CeB 229) 


The values of the coefficients, ke? COnmEne “Khas eCralt are 
listed in Table I and were obtained from the results of the 
linearization process of Appendix A. 

Speeta aememelon siould be given —o Ene coefficient 


that contains k. which is a since function of the form: 


uf 
2W Ve wiby 
Ky = [¢—E*) Sukie (ay ) J G25 oc Uy) 


aE 


Ky may be simplified by an approximation detailed in Appendix 
A of Ref. 8. The simplification is based upon the fact that 
for ahead seas, the incident frequency. ©; is; 


} V 
1l- 144 (2) 0, 


ts A= careers | aaa 22/3) 
S- 
FETE, 


For convenience in this study, the value of the ships velocity, 
Lee Hseset equal to the gravitational constant, G. This value 


of Veo mhen converted to knots, is 19.3 knots which is within 


a3 


4 
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ie.edbJe, 


> toe COLFrICiENT (xk? VALUES 


FORMULA 








CONSTANT 


24 


ARw3 VALUE 


ZAG ae Oy 


3,5336x10° 


Zoo or 


oir ne 


6.14422 


3.5052x10- 


ZOO 


(62-8319 


a) 
Ly 
(62.8313) 
Lig 
3 


8: 





the operating range of the XR-3, Continuing; 








L- Yithw. 
W, = + CR Tee ete») 
4 2 
It follows that 
I 
ie oe 20) C258 a3) 
i 2 
This approximation is valid for all Wi >2.5 radians/sec, From 
the definition of the incident wave, 
AE Ais 
a = oy (246204) 
Ky becomes 
(iP) 
Sig ZV 
Ky = A, aaa (C2 B 3 5) 
ee 
2V 
W 
Where the horizontal area of the plenum chamber is 
A, - Wohl, - C7 Bao) 
For a speed of 19.3 knots and since 
ee oe @387) 


IES) 





faeotnplifies to; 


2 
Ws 
Sin( p=) 
a 
(ats) 
ly 
Ber the XK-3, 
: 9 
A, = 200 ft 
L = 20 ft 
D 
Tinea Dee ene 
So Ce 
W 
qe = 2267, 782 saan? 
62.8319 
( ) 
W 


Ce ie) 


C2 isis SED) 


Pmeomereaduation (2.B,38), it can be seen that a null value for 


ah 


for nL = 2/2 where n is an integer. 
2. State Equation Development 


The state equations are of the form 


X = AX + Bou + GOW 


where the state variable yector 1s chosen as 


sy yt 
m 


and the output variable 


y= 0, =C X + F ow 
a 
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femwill occur for nh = Me and a peak value for k 


1 Wileee Cube 


(2.B.40) 


(2.5.41) 


C27 Bat 2) 





The states are defined as 


Xl = 8Z = (Z-Z__), Incremental Heave(Draft) CORES) 


ss 
M2 = 92 = (ZnZ_,), Incremental Heave Rate (2.B,44) 
X3 = o_ = (V_-V ), Incremental Air Mass 
5 fe outa (2.B.45) 
Mh 
where V_- = :—= (2.B.46) 
mi —. 


With control defined as a scaler, rate of incremental change 


of air mass volume, 
9 = av Ons 
m 


and where the sea input disturbance is oW. 


Beom thevsional flow graph of Fag. 2, 


= = =-J1 
xy 29 = gi t3z (2.B.49) 
a Z, 
ee A = s 7) 3s CB. 50) 
ah Z, 
: = 7 = —_ ZBreo | 
Xn a, oP Ky ( one) 


Zi 





where 


OP = COV Kank, LaaWk, -k, aZ))k, C2785?) 
Rearranging 
OP = WV Kak otk kok, Wk, kk, 02 (ee oe ee 


and grouping in order of states 


oP, = Kok3k, Wtk aks OV tk kok, OW (25825) 
The output term, y is 
Y = ae (2.B.55) 


Peeaeecuting equation (2,B,5) into equation (2.B.55) yields, 


OV -k,k,k.,k,, aW Ce 7b) 


SOUS yee Tee aan 


aXaky, dZ—k Kx 


8 3 U 


Finally, the incremental change in the volumetric air mass 


Mmemwritten as 


X3=9V_=(oUuk, , 


= -ak,k kk, d2-k kek. oy 


Zico? Ss s: m 


m ky kok 3k, Wtk, 9 du C278. 08 ) 
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Bauemions (2.B,50), (2.B,56) and (2,B.57) are the states of 


interest and are grouped below, 


fae =X? (2,B.59) 


X2 =-K pK 3k) Kg kink jk) KO X3kK KK AK, OW C2 B00) 


3 = “Ky k3k ekg X1-k3k kj X43-k) kkk, OWtK, 9 du COs Beene), 
with the output acceleration: 
y= wk kak kgAl-k jk KIX3-k) KKK, WW CZ aoero ey) 


iemomene Linear system matrices are; 


0 a 0 
A = “kyk,K) Ke 0 “K,k, Ko C22 Barts } 
“kyk3kek, 0 “k,k ko 
0 
B= 0 (2.B.64) Pe-kik,K3k, C25 si oo)) 
i 
0 
ey = “k,k,k3k, (275256) 
~k kkk, 


Ze 





Computer simulation of the linear system was accom- 
plished using the Naval Postgraduate School pregram "IODE'! 
operated from the CPCMS terminal of the IBM-360. See Appendix 
C for IODE listing. The program ocean generator for dW inputs 
a single frequency at a time using the special function fea- 
ture, Wavelength (Ly) Will change with encounter frequency 
Cw.) as shown in equation (2.B,35), It should be remembered 
that the craft speed is constant at 19,3 knots for the data 
base of this design, 

ee line Characteristic Equation 

The characteristic equation under no control condi- 

iene so: Operation, is formed from the A matrix (see Appendix 


A), 


DET(SI-A) C2mBee 7) 
S -1 0 

DEY 226 iy 8 2 S 11.248 (eBaats)) 
LEO ote 0 64.493s| 


Then the characteristic equation is; 


5° +64 .4939-+2267.782st, 5701 (C7364) 
which factors into 


(g42, SULUx107*) (9432. 2540 35,04) (S432. 25-335. 04) 


©8702 
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and 1s shown in the polenzero map of Fig. 3. Note that the 
dominant pole is a real pole close to the origin, Also note 
that the complex pair will affect early time response since 
they are relatively far removed from the origin, 

feeocot Locus plot with no contyedminput is shown 


in Fig. 4. The plot cevers a range of 


5 selon aes (22 Ba7 1) 


Table II lists the characteristic equation roots for the 


range 
o< Ke <5 0 O28. 72) 
As can be seen from equation (2.8.63), Ke 1s present 
ie 5Oth the a and a system coefficients. In order to 


om 23 
Meerorm the root locus study, the characteristic equation 


was rearranged to show the adjustable parameter Ke 


3 2 2 
Ss tk ks tka kjk kystk, (kak, koko kak otk, kjk, kok ko) tks 
(22.7 3a) 
with parameter values 
s3+(k,t10,5016)s°+2267,782s+.0928k, Gab.) 3b) 


From the table and root locus plot, it can be seen that as k, 


becomes large, the complex root pair become real with one 


Sal 





root moving toward the origin and the other migrating to 
negative infinity. Additionally, the pole near the origin 
moves very little and remains real throughout the excursion 
of K,e: it was observed that this was very similar to the 
moot focus plot in the air flow study conducted in Ref, l. 
However, the removal of the buoyancy term in this study 
caused the pole near the origin to stay real whereas the 
same pole in Ref. 1 became immaginary for large values of Ke 
ihe root locus demonstrated the high sensitivity of the com- 
plex pair to changes in air flow rate, while the pole near 
the origin exhibited a very low sensitivity to changes in 


k This observation further demonstrated that the vertical 


5° 
oscillatory characteristics of the ship depend upon changes 
tg) Kee 

As a further check on the effect of air flow on 
C.G. acceleration, an investigation of the effect of changing 
Ke on the frequency spectrum of C,G, accleration was conducted 
and the results are shown in Appendix D. From these results 
it seemed appropriate to seek a control law that would adjust 


k the air flow rate, In an effort to gain a better under- 


5 
standing of the effect of various air flow rates on the sys- 
tem, a sensitivity and stability analysis was performed in 


section IVB. The ROOTLO program used in the above study is 


shown in Appendix C, 
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Figure 3 
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iit: VERIFICATION OF THE MODEL 


A. LINEAR AND NONLINEAR SYSTEM VALIDATION 


The validation of the simplified model of a single vari- 


pleepressure lift force was conducted by comparison of its 


tPequency response characteristics with those of the model 


presented in Ref. 1. Several techniques were used to 


ae@ecomplish this, 


a; 


Single frequency constant amplitude sinusoids were 
applied to both the nonlinear system of Fig. 1 and 
also the linear system of Fig. 2. The acceleration 
output signal for the nonlinear system is shown in 
Boe. 5. Ets sintisoidal Foérm, in steady state andi- 
cates the linearity of the system, which was the case 
throughout the frequency range, 

The linear transfer function of center of gravity 
acceleration (Eq. 2.B.15) was subjected to a constant 
wave height input using a spectrum of discrete en- 
counter frequencies. The heave accelerations were 
converted to dB and plotted in a bode diagram shown 
Tee Jakes quale 

Discrete fourier transform technique, developed in 

an independent study by S, Carpenter, LRef, 9], was 
also used to test the response of the linear system 


Geis. 9)_to an irregular sea input shown in Fig. 7 
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consisting of twelve components of frequency and 
amplitude from the Pierson-Moskowitz energy density 
Spectrum given in Fig. 8 iRef. 8, 94. 

The results of procedures 2 and 3 above were compared 
with the linear system response of Ref, 1, The magnitudes 
for the center of gravity acceleration showed approximately 
20% difference between the pressure lift only, no control, 
response curve (Fig. 6) and the linear system of Ref. l 
(Fig. 10). The general shape of the curves are identical. 
Therefore, the pressure lift only model was assumed to be 
valid and useful for linear regulator control design. The 
@bserved peaks and nulls of Figs. 6, 10 are a result of the 


sinc function numerator coefficient, k as described in 


ie? 
Section III. The relationship between plenum length CL) 


and wavelength CL) expressed in k determine where, in 


oy 


the encounter frequency spectrum, the nulls and peaks occur. 
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IV, THE CONTROL LAW 


A. LINEAR REGULATOR DESIGN APPROACH 
The linear pressure lift only model is described in 


state equation form; 


xy 0 i 0 x 0 
= +. 
oS Any 0 ay. X, 0 a (4 .A.1) 
X a7 0 a5 X3 il; 
_- | 


[memerme cost function of the form 


ae 
at at 
J = v2 f (X°QX + 9U Rou) C4.A.2) 
O 
where the final value error is not included. The system may 


be designed using the linear regulator technique with the 
given quadratic cost function. An optimal control was 
found using both the transient and steady state solutions 


of the Riccati equation. The control law 
9U. z2R™~B KX (4.A,3) 
~ a ae oe 


where the K matrix wag found from the solution of the Riccati 


equation 


uy 





ek AK Oe Re BK (4,A.4) 
™ So ce 


with the appropriate choice of the Q and R coefficients, the 
feedback gains for the system were calculated. Unfortunately, 
[estinrst choice for Q and R may not yield control gains which 
are physically realizable or practical. For this reason, it 
1s appropriate to examine the effects of the control gains on 
the pole location in the s plane for the controlled system. 
Such a study was done and is presented in part B of this 
section. The choices for Q and R were derived by trial and 
error. However, the design criterion was to reduce this mag- 
nitude of the frequency response curves Pimougnie Enesstudy OF 
veetous Bode plots. To this end, many Q and R values were 
weted . 

The Riccati equation computations were done using a com- 
puter program documented in Ref. 10 and modified by B. 
DesJardins in Ref. 11. Appendix C illustrates the input/ 
output format of this program. Based upon the above design 


Gretterion, the final weighting matrix chosen was; 


‘Lo 0 0 
cen 10° 0 (4.A,5) 
0 0 5 | 
and 
R= .03 
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The optimal control gains obtained using the above Q and R 


™ ™ 


matrixes were; 


2 S 

841° -1,689x10 

Spe 1,19x10° Gore cr 
= Ik 


These gains were multiplied by their respective states and 
epmeed tO the state equations (2.8.59, 61, 62). Under state 
feedback control, the resultant output acceleration frequency 
magnitude was reduced by 13% from the peak no control magni- 
tude. Table III lists the peak lobe values of the C.G, 
Acceleration/wave height magnitudes versus frequency for the 
nonlinear and linear systems. Each system was exercised 
under control and no control conditions. The steady state 


Beausordal output for the linear model, under control and 


no control conditions is shown in Table IV Om W =H rad/sec, 
where 
= y + ee 
ete e082” 1.373 SBeesa ds, 
and 
y == a7, Cie) 


The control signal for the nonlinear system was generated by 
using the linear control low coefficient found using linear 
regulator theory. The reduction in acceleration was not 

uniform throughout the frequency range. At law frequencies, 


w <1? rad/sec, control gains caused an increase in C.G. acceleration. 
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As W approached 100 rad/sec, C,G. acceleration approached but 
never exceeded the no control condition. Effective control 


was realized in the frequency range 17<W<100 radians/sec., 


B, SENSITIVITY AND STABILITY 

Since a set of optimal gains were found for the system, 
it remained to test the system under feedback control and 
observe the sensitivity and stability for acceptability. 

The locus of roots of the characteristic equation for 
the pressure lift only system without state feedback control 
was shown in Section III, The objective of state feedback 
control is to adjust the position of the complex pair of 
roots in such a way as to reduce center of gravity accelera- 
tions. For the open control loop, the transfer function 


Peas ShOwn (Eq. 2,B.22) to be 


= Ges 1) 


Addition of state feedback will change the value of b and oc. 
Note that if the ky Sine function null effects are removed 
from the numerator, thena=c. It follows that the magnitude 
of the Bode plot will reach unity value at high frequency 
meeiratess Of changes in b and c, dn other words, the sys- 
tem has high pass frequency characteristics, Thus, the 
Besules Of increasing the yalue of c ‘will be a reduction of 
the magnitude plot in the lower range of frequencies of 


interest. Also, note that changes in the b coefficient will 


affect the resonant peaking at the corner frequency, 
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PEAK LOBE VALUES OF FREQUENCY RESPONSE, PRESSURE LIST ONLY 
LINEAR &€ NONLINEAR MODELS UNDER CONTROL AND NO CONTROL 


Lg 


CONDITIONS 
W NONLINEAR LINEAR 
NO CONTROL CONTROL NQ CONTROL . CONTROL 

6 ile eso Zomeo ISAS Deo 
y 1s as Zoe Itch. 254 5 
8 horas Zia Seer 21.4 
9 9.7 1 a 9,8 ey 
14 Lae 42.2 42.3 45.7 
dats 48,2 48,2 48.4 5029 
16 49.7 yore 7 49,9 Sse 
aay u5.9 45,9 46.1 46.2 
24 69.8 62.9 TORS 63.4 
25 Wicca 69,7 79.0 0 ea 
26 (seal 70 79.9 TONS: 
27 liga a4 fife xis) 63.9 
28 

34 90.8 Vi eal 91.8 (Sas 
50 LEOORS S53 INI. S723 
36 W244 86.8 MO) 258 87.9 
37 Se 78.6 92.6 79.6 
38 

44 99°5 86.8 TAG = oe By ee 
45 Lae Oneme, 96.4 eS 0) Sis) |e, 
46 MEO 0 96,8 cers 98.0 
47 99.7 88,0 Ielsl 88.9 
48 

54 97.1 88.5 97,5 88.5 
3:5 Gy ef 98,5 Oe 0 Shoe i 
Se ROS. L 99.2 IOS ees OG, 0 
5) 9) 98.2 S045 99.0 Spee 








TABEE £it 
Ccontinued) 


NONLINEAR 
NO CONTROL CONTROL 
88.7 83,8 
99.3 gu 1 
100.4 95,4 
91.9 87.6 
79,0 76.5 
89.5 86.9 
91.3 88,8 
guy 82,2 
69,9 68.8 
80.2 78.9 
82.6 81.4 
aes 76.1 
61.9 61.4 
Pil, @ fale 
75.0 74.5 
70.9 70.5 


u9 
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ee 
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TABLE IV 


SINUSOIDAL LINEAR MODEL QUTPUT WITH/WITHOUT PEEDBACK CONTROL 
AT Ww F45,0 Rad/Sec. 


TIME CONTROL CONTROL EFFORT NO CONTROL 
i a : 
0.780 54,38 Hees 61.02 
0.785 34.33 2,64 38.03 
0,790 12,50 0.86 13.08 
0.795 -9.96 a. 95 pilose 
0.800 =e Be Be <6), 89 
0.805 B52 .22 Bineas -60.56 
0.810 -69,85 ~5.76 ~80.49 
0.815 ~83,86 BY; See 
0.820 AGIs 5 3G BY. -107.08 
0.825 ~98 44 7 ee 2 pee 
0.830 So) me 29) 3. Sak. 8 
0.835 =S PAu ag alae -105.62 
0.840 =Sesi0 -~6.64 =o oy 
0.845 -~68.60 =5 45 =97 281 
0.850 Geran ES ecor -56.85 
0.855 2310427 EY: S255 6 
0.860 act Hi) SS Bou 
0.865 rue up, 2 17.40 
0.870 35.95 3.04 42.10 
0.875 55.82 u.63 64.64 
0.880 72.82 5.99 83.85 
0.885 86.06 TOs 98.73 
0.890 94.87 AUD 108.53 
0.895 98.79 8.00 Oe, 72 
0.900 97,61 7.88 111.11 
0.905 ae ae 103.79 
0.910 80,50 6.43 27. Lal 
0.915 65.44 5.19 73.74 
0.920 u7.01 3.68 52.56 
0.925 26.49 1.98 28.68 
0.930 3.94 Dome goal 
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lia The Characteristic Equation For Feedback Control 
— es 


The signal flow graph of Fig. 2 was modified to show 
Feedback gains and is shown in Fig. 11. Using Masons 
gain rule, the characteristic equation was derived from the 
S.F.G, (Eq. 4,B.10). The resulting equation was compared 
to characteristic equation derived from the modified system 
matrix with feedback of Appendix D. It was noted that they 


differ by the terms 


and 45 3435 


rele) 
These terms very nearly cancel and are insignificant (a5) 
Ag 3745 3437,7-970D when compared with the other larger terms 
in the characteristic equation. Therefore, their omission 


maemmeiie signal flow graph of Fig. 11 is valid: The deri- 


weoaerom £LOllows: 
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=— i +L, + 
A LL, L +L, Gi dlie: * Ly L. Cur BT) 
g a a a = g =a 
ec ee 23 23 3 21 
S 9 S S 
(4,B.8) 
dh 3 Z. 
ed ah a 
A 5 \s Z4S + A448 + A578 + 253858 25384 
2 Bay (4,B.9) 
s?t(a +g ys? Geel. Oneal oe tea oe saeco) 
a5) 6 Pies EE 2123 Vig ag 
(iB Lo) 


weeeOOc Locus 

The following root locus plots are a study of the 
effects of changing the three different control feedback gains. 
Many weighting combinations were tried and it was discovered 
that the system was most sensitive to changes in air flow rates. 
This observation was in agreement with the analysis of Section 
II that the system was most sensitive to kee Each of the 
following root locus plots and associated table of values was 
chosen as representative of the alternatives tried, The 
characteristic equation used is shown as Eq, (4.B.10). The 
Naval Postgraduate School program for the IBM 360/67 computer 
was used to generate the tables and plots. The program 


listing is shown in Appendix C, 


Ss) 





In the first study, the constant gains used were 


= U7 5), 0 =a 


S11 » 812 


and 8,3 was varied over three decades starting from 843 > alee 
The resultant plot and table of values are shown in Fig. 12 
and Table IV respectively. For ease of reference, the sys: 
tem poles were given letter designations for the characteris- 


tic equation 


(StA)(S+B)(S+#C) = 0 (4.B.11) 
where 

pet 0? 25-425. 00) 

B = (+32.25+535.04) 

c= (2.5x107'+40) 


As 8,3 nears -20, pole a rapidly enters the real axis and 
migrates towards negative infinity and pole B moves upward 
across the real axis, Pole C moves in a negative direction 
Mie aL 8137726 and then maves below the. real axis tracking the 
mirror image path of pole B. It was seen from this plot 
that as B,3 nears -20 that the complex pair undergoes the 
greatest movement. In this study there was no point at 
which instability occurred, When the gains were tested on 
the system, there was no noticable decrease in vertical 
acceleration. 

The next phase of this investigation uses constant 


gains of 
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while B13 varied over two decades from -.1. These gains 


were calculated from the Riccati with weighting matrix; 


106 0 0 
oe 0 a 0 (4.B.12) 
0 0 5 
and 
i = 03 Gu pels) 


The resulting plot and table of pole values are shown in 
Fig. 13 and Table V respectively. As B33 varies, the poles 
follow similar paths as in the previous study. The greatest 
degree of pole movement, or sensitivity, occurs at 8, 37-16.0. 
This set of gains produced instability when 8,37 71.08. 

In the final study, the constant gains were 
= -1688.095, 


11S eG. = -18.306 
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Again, B13 was made to vary over three decades starting with 
813 7 -~.l1. The resultant plot and table of values are shown 
in Fig. 14 and Table VI respectively. These gains were cal- 
culated using the Q and & values shown in (4,A.5). Again, 

the poles follow similar paths but pass closer to the origin 


than in the previous two studies, causing a greater net 


decrease in the systems oscillatory tendencies. As in the 
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previous study, some values of 213 caused instability. 
mewtacularly, for B437799, pole C was in the right half 
plane, The greatest sensitivity was experienced for values 
of B13 near g,3 = ~30. Note that the value of 813 chosen 
via Riccati will cause stable operation and is also removed 
from the value of greatest sensitivity. 
2. kouth Criterion 

The Routh criterion was applied in order to confirm 
the range of feedback gains for stable operation obtained 
from the root locus study. The Routh array is set up below. 
Mmiemenee., fiven in Eq. 4.B.17, is used with the understanding 


that all coefficients are absolute values, 


. oO eB 29 
433" 83 493 21 * 254 837 
gan 723 212 ~ Ss 0 

S538] ° 72153 v 
Diemweonditions for stability are: 
IL} ZB, 7 7433 (4,B.14) 

since 4,3 - 64.493, then 

Za -64.493 for stable operations. 
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Goals oe Roe, 


(4.B.16) 





If g, > ~64 the g, > -201,72. 


The boundary limits of g,, and therefore g,, are 
very sensitive to Ba: 

For the specifie control chosen in Section IVA, the 
Routh array was used to determine the stability limit for 
the gains, 


meom Eqt,B.14) 


213 33 
therefore Tite) Sy eae ele Se Cael: 
Baomesg( 4.3.15) 

ia me $13 
23 


therefore foie oma soo o.oo. 
Emeom bq(4.B.16) 


( ) 


Ss Shae | Seis 
- a 
é 53643378) 3 
therefore tee > =72/3295 


It was seen from this study that the gains were well 
within the limits of stability. Lower gains were realitively 


more stable. 
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V, CONCLUSION 


It has been shown that the linear regulator design can 
be applied to the problem of vertical plane motion for the 
captured air bubble surface effect ship. The linear pressure 
lift only model, developed from the simplified equations of 
motion, was shown to be a valid model which could be used for 
control analysis and design. An optimal control was found, 
based upon a quadratic cost function, using empirically de- 
rived weighting factors, The pressure lift only linear 
system showed a reduction in center of gravity acceleration 
for 17 <w <100 radians/sec. A C.G, acceleration reduction 
was observed in both the linear and nonlinear systems, The 
shape of the frequency response magnitude curves were identi- 
cal. A 13% maximum reduction was observed at Deo rad/sec. 
where peak acceleration occurs. 

The control gains derived, as shown in Section IV, are 
Within the airflow rates currently available for the XR-3. 
The Riccati equation computer program on Ref. 10 proved to 
be a most useful tool, which facilitated experimentation of 
a large number of weighting coefficients, making an optimal 
control design more feasible, 

Finally, the analysis shows that the system is stable 
over the range for the control gains selected. It was seen 
from the root locus study that the limits of stability narrowed, 


Mere Sensitivity to air flow rates increased and C.G. 


GS 





acceleration decreased with higher control gains. It was 
seen that for these higher gains the C,G, acceleration was 
not reduced and the control, 94? increased fon we <l7 rad/sec. 
This was due to the fact that the buoyancy effects are pre- 
dominant in that frequency range and the buoyancy terms are 
at a constant steady state value in this model. 

Though an optimal control was found based on a particu- 
lar quadratic cost function, many other possibilities exist. 
Different defined weighting factors can be experimented with, 
and some may result in improved C.G, reduction. The limiting 
factor in such a design is equipment limitation. The choice 


of weighting factors, and ultimately, the control, are even- 


tually determined by that fact. 
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YI. RECOMMENDATIONS FOR PURTHER STUDY 


The following list of recommendations is presented for 


consideration and further study. 


A. A more realistic system presentation could be accom- 
plished 1f pitch coupling was combined into the 
linear regulator study. This would facilitate the 
speed analysis study for the avoidance of wet decks 
and undesirable pitch accelerations. 

B. The passive control technique discussed by Boggio 
in Ref. 6 should be examined further. The introduc- 
tion of a flexible membrane into the plenum chamber 
may reduce vertical accelerations. This technique 
has been examined in model towing tank studies. 

C. Since the mass of plenum air X3 cannot be directly 
measured, an observer might be used in conjunction 
with state X3 or perhaps remodel using low order 
model analysis. Perhaps consideration should be 
given to a design using only X1 and X2 as states to 
be used for feedback, 

D, The linear regulator technique design approach should 
be compared in greater detail to the filter design 
technique proposed by Rohr Marine Inc, in Ref. 3. A 
fixed air flow rate, based upon a filtered observed 


sea state, could be a more realizable design in terms 
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of the physical equipment limitations. The study 
of fixed air flow, of Appendix E, shows that the 
friter design approach may offer greater control 
over the entire encounter frequency range. (see 


Dabfege Eal) 
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APPENDIX A 


The simplified non linear system of Ref. 1 is shown in 
Fig, Al. This system, with the buoyancy term removed, served 
as the basis for study in this thesis, The linearization 
and resultant linear system matrix is shown for the simpli- 
fied non linear system, The pressure lift only model lineari- 


Pouerom 1S also shown. 


mojeeaetZATTON OF THE SIMPLIFIED MODEL 
The equations of motion are linearized using Taylor 
series expansion around the steady state operating point. 
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Writing the above equations in state equation form: 
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Uwe characteristic equation was found from the system matrix, 


DETCsInA) = 0 


a3 





=o De 
5 
Soe Gos + 283s nh 739 = 0 


piem factors into; 
(24137) (st 30. 624724, 8) (st30,82-934,8) 


TAYLOR SERIES EXPANSION OF PRESSURE LIFT ONLY VERTICAL 
EQUATIONS OF MOTION 
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then incremental acceleration due to pressure nS 
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Defining the air mass volume: 
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Then the incremental volumetric change in air mass is; 
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and the incremental force due to pressure is; 


Incremental air flow into the plenum 
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Incremental air flow out of the plenum 
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Incremental gauge pressure (plenum) 
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eeyeremental Pressure (plunum) 
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Plenum Pressure 
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A series of system coefficients, k 2» were defined for ease 


of writing and are used throughout the thesis. They are, 
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APPENDIX B 


A brief step weight removal study was performed and 
compared ta Ref. 1 as an additional check on system validity. 
The IODE remote terminal CP-CMS program on the IBM/360 was 
used in this study. A sample listing is shown in Appendix 


C. 
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eee ECTPICATIONasa—— 


VARIABLES €& INITIAL CONDITIONS; 


X1 = 0.500000 
X2 = 0,0 
X3 = -44.5300 
i= 0,0 


SPECIAL FUNCTIONS: 
Y = 2030.26*X1 10.09*x3 


DERIVATIVES: 

mexi)/D(T ) = = 
X2 

Meee /D(T ) = = 
-2182.68*X1-10.76*X3 

D(X3 /D(T 0 = = 
-12504.36*X1-62.0*X3 


CULeUTS : 
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SO meAeCULATION WHEN T .GE. 0.600000 
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APPENDIX D 


The characteristic equation for the pressure lift only 
model under control conditions was derived from the system 
matrix with the feedback control gains added to the state 


Xa The result was compared with the characteristic equation 


Srecection IV.B.1. 


ol: 





The system matrix is 


0 1 0 
a= any 0 a5 3 

a 0 a 

were a6 


and with feedback gains; 


0 1 0 


ot 23 


+ + 
ee eo B12 4337813 


The characteristic equation with gains becomes 


3 2 
= + =e 
S$ ~(43348)3)8 (4538) %a51)8 
“ele een 
7514337491813 


It should be noted that the above equation is identical to 


that derived in Ref, 1 except that the gains are shown in 


yee.) 
any te Ao ane 
An 4 = all .248 
Az, =I5002 28 
a> = ~6§64.493 


g2 





From the Riccati equation, it was seen that the gains had the 


following sign values; 


Tio ole 2 


substitution of the appropriate signs yields 


3 2 
S +(a337813)8 +(ay38)9ta51)5 


749743 37991813 


PG eset Day eal 


Substituting in the system value parameters, the character- 


istic equation becomes: 
3 2 
S +(64.4934tg),)s #(2267.782+11, 248g, ,)s-1l.248g,, 


ee ieee 5710 
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Ar PEND e 


EEXED ATR FLOW RATE STUDY 


A brief study was conducted to examine the effects of Ke 
on the overall frequency response and compared with the op- 
timal control design results, Earlier work indicated that 
the system response was sensitive air flow rate. Therefore 
three values of Ke were tested on the pressure lift only 


non linear system with the resultant response curves plotted. 


Values for k. are: 


6 
ke = 6.144 , Fig. El 
ieee) 12.0 , Fig. E2 
i= 29), 0 m Fig. E3 


Note that the high acceleration midband frequencies are 
attenuated as Ke increases. The response curve flattens 

out and midband frequencies are reduced by 20%. This result 
would tend to confirm the contention, of Ref. 3, that fil- 
tered fixed air flow rates method of control may have 
advantages. The time response for sinusoidal input remains 
linear throughout the range of interest, The time versus 
acceleration plot shown is a time history for an air flow 


mace Of Qr0 = 60 at an encounter frequency of Waa rad/sec, 


Appendix C presents the DSL Bode program listing. Finally, 


gu 





the system peak frequency magnitude plot for the conditions 
of no control, state yariable feedback control and fixed air 


Etow Date control is shown in Fig, E8 for Q70780- 
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